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SUMMARY 

Sounding rocket Experiment 74-18 was designed to demonstrate the 
effects of the Black Brandt rocket acceleration levels (during the low-g coast 
phase of its flight) on the motion of a liquid metal system. The experiment 
Principal Investigator was Mr. Charles S. Schafer of NASA-MSFC. This 
Lockheed -Huntsville study was performed primarily to assist the Principal 
Investigator in preflight design. Some post-flight analyses were also conducted. 

Preflight studies consisted of heat transfer analysis and convection 
sensitivity and convection modeling analyses. These studies aided in the: (1) 
final selection of fluid materials (indium-lead melts rather than paraffins); 

(2) design and timing of heater and quench system; and (3) preflight predictions 
of the degree of lead penetration into the pure indium segment of the fluid. 
Postflight studies involved: (1) updating the convection sensitivity calculations 
by utilizing actual flight gravity levels; and (2) modeling the mixing in the 
flight samples. 

Preflight convection sensitivity calculations were made for the cylinder 
perpendicular and parallel with the acceleration of gravity. The values for 
the Marangoni, Grashof, Praudtl, Rayleigh, Nusselt, solutal Rayleigh and 
Schmidt numbers, momentum diffusion time, and thermal conduction time 
are calculated for the melted indium. These sensitivity calculations indi- 
cate that: 

• Expected convective heat losses, from the mixing cylinders to the 
ambient air, during rocket tests can be simulated by ground tests 
utilizing ambient air at 0.1 atm. 

• Thermogravitational (including g-jitter) and thermocapillary con- 
vection will not induce mixing during the rocket tests. 

• Solutogravitationail and possibly solutocapillary convection will 
induce mixing during rocket tests. 
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In addition to sensitivity analyses, the heat transfer in the mixing experi- 
mental apparatus was mathematically modeled via a lumped parameter model 
(Lockheed-Huntsville Thermal Analyzer Computer Program). Graphs of 
the heat-up and quench were drawn which show the temperature history of 
the indium/indium-lead mixture. Based on these results, recommendations 
were made to move the quench time from 310 to 295 seconds (laboratory time 
point). The temperature gradients calculated were also used in the convection 
sensitivity studies. 

The postflight sensitivity analyses confirmed rocket results because 
they indicated that Flight Sample No. 6 was most susceptible to convection; 
while the other two samples(Nos. 1 and 2) were subject to only marginal or 
no convective mixing. The analyses indicated that only solutal (concentration 
gradient) induced buoyancy and surface tension forces can induce stirring in 
the flight samples. The postflight modeling results, however, did not show 
significant mixing in any of the flight cases. These models only considered 
buoyancy forces. Thus, it is concluded that the large amount of mixing seen 
in Flight Sample No, 6 is the result of solutal surface tension forces or the result 
of unexpected solutal buoyancy forces. The latter is most probable and can be 
explained if the flight acceleration levels were higher than those used in this 
study and/or if a more rigorous (three-dimensional) computer model is 
utilized. It is recommended that further analysis of acceleration levels and 
that several additional cases of the present computer model be conducted to 
clarify the rocket flight results. 
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Section 1 

PREFLIGHT ANALYSES 

This section presents the results of preflight thermal analysis and 
convection prediction portions of a theoretical investigation into the Indium 
and Indium-Lead Mixing Experiment for the Space Processing Rocket Program 
of NASA-MSFC (designated as Experiment 74-18). The phenomena being 
investigated are those which could induce flow by causes other than residual 
acceleration through the non- spin platform of the rocket system. In this 
first phase investigations were specifically directed to supporting the design, 
fabrication and flight sequence timing of the Mixing Experiment 74-18. Other 
theoretical phases of this program include convection calculations from various 
causes such as thermal gradients, surface tension, bubble convection, phase 
change and g-jitter. The results presented in this section include temperature 
histories of the metals during heat-up and quench and convection sensitivity 
calculations. 

1.1 DESCRIPTION OF THE APPARATUS 

The experimental system is a set of mutually perpendicular cylindrical 
samples which were melted, allowed to remain molten for about four minutes, 
then resolidified, all while in the low-g portion of the rocket flight. The 
samples were composed of two cylindrical sections joined together; one 
section was pure indium, the other was 80% indium and 20% lead (by weight). 
This gives a density difference of about 7% between the two liquids. 

The mixing apparatus is basically an aluminum cylinder 0.614 cm i.d. 
and 3.68 cm lonf which is closed and sealed at both ends by aluminum caps. 
One-half of the free length of 2.5 cm is filled with indium metal except for an 


Engineering drawings of the mixing apparatus are shown in Figs. A- 1,-2 and -3. 
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air space at the end for expansion. The other half is filled with an alloy of 
indium 20% lead except for the expansion air space. The inner end surfaces 
of the metal and the alloy touch at the midpoint along the length of the cylinder 
and are held in position by "smoked" aluminum rings at each end of the 
cylinder. 

The indium and indium-lead are heated past the melting point by a ni- 
chrome heater wire wound around the cylinder. This heater wire has a stain- 
less steel sheath and is electrically insulated from the sheath by magnesium 
oxide powder. There are 29j turns wound around the cylinder and the heater 
wire then goes directly to a terminal on the mounting bulkhead. 

The cylinder is held in position and supported by a thin stainless steel 
bracket which is welded to the heater wire sheath and attached to the mounting 
bulkhead. Also mounted to this bulkhead are two other (identical) cylinders 
whose axes are perpendicular and are also mutually perpendicular to the 
first cylinder's axis. These two cylinders are supported by a single mount- 
ing bracket, while the other cylinder is mounted by a single bracket. These 
three cylinders, heaters, mounting brackets and bulkhead are contained inside 
a cylindrical chamber with a void volume of approximately 571.8 cc. 

In order to resolidify the metal and alloy very quickly a supply of water 
is carried in an elastic bladder attached at one end of the cylindrical chamber. 
The contained water is released at 295 sec (laboratory time) after the heater 
is turned on, and is rapidly expelled into the chamber by water pressure 
which is created by the elasticity of the bladder. There is an excess of water 
in the bladder which ensures that a large percentage of the chamber will be 
filled. An air relief tube vents the air behind the bladder as the chamber is 
being filled. The water which is filling the chamber, is in motion and exposes 
the heaters to fresh cold water almost continuously. This will quench the 
molten metal and alloy before the rocket payload experience high-g forces 
upon reentry. 

There are two time frames in which this experiment can be investigated; 
one is the rocket time frame; the other is the laboratory time frame. These 
are compared in the following table. 
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Event After Rocket 
Firing 

Rocket Time 
Frame (sec) 

Laboratory Time 
Frame (sec) 

Heat Turned On 

45 

0 

Heater Turned Off 

225 

180 

Water Quench Initiated 

340 

295 

High-g Experienced 

405 

360 


All times given in this report will be in the laboratory time frame unless 
otherwise specified. 

1.2 THEORETICAL THERMAL ANALYSIS 

In order to investigate the heat transmission from the heater to the 
indium metal and the indium -lead alloy a mathematical model of the mixing 
apparatus had to be made. The partial differential equations to be solved 
which would describe the transmission of heat by conduction through the 
apparatus, are sufficiently complicated that a discrete model and computer 
solution is needed. The best type of discrete model is the lumped parameter 
model. From the engineering drawings a lumped parameter model was devel- 
oped which considered the geometry, materials and heater of the mixing appa- 
ratus. The math model was made with the least amount of deviation from the 
actual apparatus. The correct dimensions, density, thermal conductivity, and 
specific heat for each material were used. In the case of the indium and indium- 
lead and the surrounding air, each were dimensionally subdivided in order to 
obtain better thermal gradients and conduction patterns. The lumped parameter 
modeling method assumes that the model is divided into cells so that the prop- 
erties of a particular cell may be stated as existing at the center of the cell 
or at its surfaces. This method leads itself to the solution of complicated 
problems that involve partial differential equations which describe the behavior 
of physical parameters. Thermal conduction is one such problem and the 
temperature is a parameter of interest. 

The model has the indium solid cylinder divided into six pieces; two 
radial sections by three longitudinal sections. The indium-lead alloy is also 
subdivided in a like manner. An air space has been left at the outer ends of 
the solid cylinders for expansion. The aluminum end caps have each been 
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subdivided into two longitudinal pieces. The aluminum cylinder was not 
subdivided and is considered to have exactly the same dimensions as shown 
in the engineering drawings. 

The heater wire which is wound around the cylinder Z9j turns has been 
modeled as a set of concentric thin cylinders. The heater wire is nichrome 
with a stainless steel sheath and a magnesium oxide electrical insulation 
powder packed in between the nichrome wire and the sheath. Since the stain- 
less steel sheath is very thin, the heat is conducted through it fast, so it can be 
represented as a cell of zero thickness. In order to account for the distance the 
thickness of the sheath is added to the thickness of the magnesium oxide. This 
simplification of the heater model was checked by leaving the sheath in the cal- 
culations and comparing the resulting temperatures, which were almost the 
same in both cases. So, the heater was modeled as three concentric cylinders; 
the inner one of MgO; the center one of nichrome; the outer one of MgO. 

The air was subdivided into a number of concentric cylindrical rings 
around the heater and at the ends of the cylinder and heater. The air ex- 
tended 2.5 cm in an axial direction from the end of the cylinder and heater, 
and it extends a total of 3 cm radially from the centerline of the cylinders. A 
"cold wall" (22°C) heat sink was placed at the outer edge of the air and com- 
pletely surrounding all the air. 

During the testing portion of the lumped parameter model, a number of 
different computer trials were made using various model configurations. In 
one such trial the holding bracket was included; in another the radiation was 
considered. In still a different trial am. effective thermal conductivity for 
air was used to consider convection at Ig. Since one of the ground tests 
was run in a vacuum chamber at 0.10 atmosphere of pressure, a computer 
trial was used to approximate the convective-conductive air at this pressure. 

In this case an effective thermal conductivity was calculated and used which 
considered convection at 0.10 atmosphere and the thermal conduction in air. 

To effect the water quench the computer program was stopped at the 
295 second step, and values for water were substituted for the air. With the 
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air replaced by water which was a heat sink (because the temperature wouH 
only rise 3.8°C) the model was run on the computer until the 360 secorn tep. 

In the lumped parameter model the temperature of a particular ceil is 
taken at the center of that cell. Also considered at the center of the cell is 
the thermal capacity of the cell's volume. To propagate the heat by conduction 
each cell is connected by a heat flow resistance. This resistance is inversely 
proportional to the thermal conductivity of the cell and inversely proportional 
to the bounding area of the cell through which the heat must be conducted. This 
resistance is directly proportional to the path length which is the distance 
between centers of the two connecting cells. Radiation is directly proportional 
to the fourth power of the temperature differences. Thus, the parameters of 
the lumped model are defined at the center of the cell, at the surface of the cell, 
and by cell interconnecting parameters. 

Physical problems solved through a lumped parameter model yield real- 
istic solutions that match closely to experimental results. The lumped param- 
eter model is illustrated as Fig. 1. The thermal properties of the materials 
in the mixing apparatus are given in Tables 1 and 2 (utilizing Refs. 2, 3, 4, 5). 

1.3 THERMAL ANALYZER COMPUTER PROGRAM 

Prior to this project Lockheed 'Huntsville developed a computer program 
which was used to solve heat transfer systems. The Thermal Analyzer com- 
puter program utilizes a finite difference technique to solve the general heat 
balance equation for systems in zero or low gravity. The lumped parameter 
model thereby lends itself to easy solution by the Thermal Analyzer program. 


The general heat balance equation that this computer program solves for 
special cases is (Ref. 1); 


'^i,e+ AG 




R.. ^i,G + AG 
U 


AG 

C. 
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Fig. 1 - Illustration of the Lumped Parameter Model 


PAGE IS 
POOR QUALITY 


LOCKHEED- HUNTSVILLE NESEANCH i ENGINEERING CENTER 


i 




LMSC-HREC TR D^9<)H46 


Table 1 

THERMODYNAMIC PROPERTIES* OF THE MATERIAES 
IN THE MIXINr APPARATUS 


Material 

P 

k 

C 

P 

In(156.4°C M.P.) 

7.25 

0.06 

0.058 

In-Pb(20%) (160°C M.P.) 

8.07 

0.064 

0.053 

In (liquid) 

7.02 

0.10 

0.0654 

In-Pb (20%) (liquid) 

7.82 

0.088 

0.0596 

A1 6061-T6 

2.80 

1 

I 0.35 

1 ! 

0.211 

Quartz 

2.65 

0.023 

[ 0.0165 

Cu 

8.94 

0.948 

0.0922 

Nichrome 

8.247 

0.032 1 

i 

0.107 

MgO 

3.58 j 

0.0002 

0.24 

Stainless Steel 347 

7.85 

0.088 

0.125 

Stamless Steel 304 

7.85 

0.088 

1 0.125 

1 

Air 

0.0012 

62.2 X 10"^ 

0.238 

Water 

1.0 

0.0015 

0.998 

Eicosane (36.8°C) 

0.841 

35.6 X 10"^ 

0.530 

Octacosane (61.4°C) 

0.912 

35.6 X 10"^ 

0.511 

Eicosane (liquid) 

0.744 


I 

i 

i 

Octacosane (liquid) 

0.756 

1 

1 
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Table 2 

PHYSICAL PROPERTIES OF UQUID INDIUM 
AND INDIUM- LEAD 


Property vomits) 

Symbol 

Value* 
for In 

Value* 
for In-Pb 

Dynamic Viscosity (poise) 
(gm/cm/ sec) 

M 

0.012 


Kinematic Viscosity (stokes) 
(cm /sec) 


1.7x10'^ 

{ 

Surface Tension (dyne/cm) 

o 

340.0 


Change of o with Temperature 

. .Q 

(dyne/cm/ C) 

do 

dT 

- 0.08 


Volumetric Expansion Coefficient 
(cm /cm / C) 

P 

LISkIO"** 


Thermal Diffusivity (cm /sec) 

a 

0.218 


Impurity Concentration (mole/mole) 

C 


0.20 - 0.23 

Change of Density with Concentration 

Pc p \9C/^ 

0.497 


2 

Mass Diffusivity (cm /sec) 


1-3x10"^ 


Heat of Fusion (cal/gm) 


6.8 

6.6 

Radius of In (or In-Pb) (cm) 

r(H) 

0.307 

0.307 

Length of In and In-Pb (cm) 

U(H) 

2.35 

1 

1 

1 

Mass (gm) 

pv 

2.52 

2.80 

Gravity at 10"^g (cm/sec^) 

g 


9.8 X 10"^ 


All values are given in CGS units. 
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whore 

th 

= tlie temperature for the i ” eell 
0 = the time variable 


A0 = 

R.. = 
ij 


time increment 

resistance to thermal conduction between the i^^ cell and 
the cell. 


= average heat rate over the time increment A0 (cal/sec) 
(heat added from the heater or taken up by a heat sink) 

= heat cajjacitance of the i^^ cell 


where 


and 


R.. 

ij 


L 

kA 


C. = pVC 


L 

k 

A 

P 

C 

,.P 


= effective condo ~tion path length between cells (cm) 

= thermal conductivity of the cell's material 

= effective cross-sectional area for heat transfer between two cells 
= density of the cell's material 
= specific heat of the cell's material 

3 

= volume of the cell (cm ) 


Radiative heat transfer computed from equation: 


Q 


rad 




where 



2 o 4 

Stefan -Boltzmann constant (cal/s ec-cm -( K) ) 

2 

effective radiating area (cm ) 
effective view and emissivity factors 
absolute lemperature ot me i ceii ( x\> 
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The radiative resistance is calculated by 


A c^.(T^ + T^)(T. + T.) 
r 1 J ' ' 1 J 

which is used in the resistive parts of the general heat transfer equation. 

The Thermal Analyzer program computes new temperatures for each 
cell in time step increments. An iterative method is used to achieve accurate 
results in solving the heat transfer equation for each cell. Thus, very accu- 
rate temperatures can be calculated for each cell of the lumped parameter 
model at certain time increments. The time increments are computed by 
the program and are determined by the stability of the solution. In this par- 
ticular case of the mixing apparatus, the time step was 0.0015 sec, or 667 
steps to equal 1 second of the laboratory time. 

1.4 THERMAL RESULTS 

Entering the lumped parameter model into the Thermal Analyzer pro- 
gram yielr's a time -temperature trace of the mixing apparatus. The program 
calculates and prints out the temperature of each and every cell, from 0 sec 
to 360 sec in increments of 5 sec. 

The heat rate put out by the heater was hand-calculated by the equation 

° = OfR 

where the applied voltage on the heater was Z8 volts and the heater resistance 
for one heater is 46.5 ohms. This yields 4.03 calories per second heat input 
from the heater wire. This quantity of heat was input in the nichrome wire 
of the lumped parameter model. 
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The phase change at the melt time was hand calculated so that fewer 
computer runs could be made. The calories required to melt the indium 
metal and the indium-lead alloy is: 

qi„ = qi„.pb= (P''Hf)i„-Pb 

where 

p = density 
V = volume 

= latent heat of fusion 

The subscript In and In-Pb designate quantities used in the paranthesis are 
for indium and indium-lead respectively. Assume that one-half of the heat 
input goes to the air and heat sink (which is close to that value computed by 
the computer program), and that all influx of heat at the time of melt goes to 
melting the metal. The time for melting is calculated from 



where 

Q = heat rate input from the heater wire. 

These calculations yield a melt time for the indium of 8.49 sec and a melt 
time for the indium-lead of 9.17 seconds. Since there is approximately 3.6°C 
between the melt points and thus only 2 seconds of duration occur after the 
indium melts, the two melt times can be additive. The time element between 
melts is also shortened by heat input to both so that the indium-lead alloy will 
probably begin melting before the pure indium has completely melted. These 
effects would tend to cancel thus yielding a total melt time of approximately 
17 seconds. Only the melt time for the indium is plotted in the following 
graphs. 
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These graphs, shown in Fig, 2 ., a.ro the temperature versus time plots for 
jidium only. The graphs arc for the cell which is located at the outer edge 
near the air gap and at the center of the cylindrical cell. The hand calculated 
melt time for indium is added to the computer calculated temperature history 
to yield a complete curve. The curve for the indium lead is almost identical 
except for the melt temperature. 

i.5 TEMPERATURE VS TIME GRAPHS 

Three different cases were investigated: (1) low loss case; (2) typical 
case; and (3) highest heat leak possible case. These different cases were con- 
sidered because cases (1) and ( 3 ) should bracket the actual flight experiment, and 
the typical case lies between them. The low loss case considers the smallest 
heat leak that the experiment will probably have. For this minimal heat loss 
case only the bulkhead and the container are considered but they are considered 
not to gain heat. Thus, a "cold wall" heat sink at 22°C is considered to sur- 
round the heater and cylinder at the outer edge of the air. In this case the only 
heat lost is that which is conducted through the non -convective air to the wall. 

The typical case considers the "cold wall" along with a holding bracket 
and radiation to the wall of the container. The bracket in the typical case was 
modeled as a combination of the two brackets used to hold the three different 
mixing cylinders and heaters. One of the actual brackets is longer than the 
ochtir and it holds two of the cylinders while the shorter one holds only one 
cylinder. I diation losses are also considered from the exterior of the heater 
and cylii ’er to the bulkhead and the container walls which are 22*^C. The heat 

'if, 

loss values are: 


• "Cold Wall" - 2% 

• Through the Bracket ~ 3% 

• By Radiation to the "Cold Wall" ~ 4%. 

The highest heat leak case considered the "cold wall," bracket and radia- 
tion again, but the bracket in this case was shortened and mounted on a heat 


Fraction of convective heat loss to surrounding air. 
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Time the heater is turned "on" 
Time the heater is turned "off 
Time the quench is started 


Time that high- 
Temperature of 
Temperature of 
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sink bulkhead. This bracket was modeled after the shortest of the two brackets 
which only holds one heater and cylinder. The bulkhead to which it is attached, 
was considered as a heat sink at 22°C. The heat loss rate through this bracket 
was approximately at a 10% rate, which is 2.5 times more than that lost from 
the typical case bracket. Again radiation and the "cold wall" contributed 4 and 
2%, respectively, to the heat loss rate. 

The water quench was modeled by substituting water for the air. This 
was accomplished by substituting water specific heat, conductivity and density 
for those same properties for air. The temperature of the water was set and 
held at 22°C which made the water a heat sink. The initial temperature of 
each cell, except the cells occupied by water, was taken from the heat-up, 
computer run at 295 sec of lab time. There was a quench computation made 
for each of the three cases. The duration of the phase change during the 
quench was calculated from the slope, which yielded 3.8 sec. 

One further case that was computed was that of convective air at 0.10 
atmospheres of pressure. This case was computed in order to match the 
conditions of a ground test which was made in a vacuum of 0.10 atmospher .s. 
The Nusselt number for 1 g and air density of 0.1 p was calculated and yielded 
3.88 which is 31.6% of that of air at one atmosphere and 1 g. Since the Nusselt 
number is the ratio of the convective heat transfer to the conductive heat trans- 
fer, the thermal conductivity of the air was modified by a multiplicative factor 
of (1 3.88). The air then has it's usual thermal conductivity plus the convec- 

tion due to 1 g and a pressure 0.10 atmosphere. The temperature-time graph 
(Fig. 3) of this case matched very well the data obtained in the vacuum ground 
test. The test data are not yet available in graph form, thus the data are not 
shown in Fig. 3. This vacuum ground test was conducted to simulate convective 
heat losses (from the cylinders to the ambient air) during rocket testing. 

1 .6 THERMAL GRADIENTS 

The graphs shown in Figs. 2 and 3 are for one typical cell of indium; how- 
ever, each cell of the indium and the indium-lead is very near this temperature 
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E 3 Temperature of Indium melt point 
F 3 Temperature of Indium-20% Lead melt point 
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Fig. 3 - Time-Dependent Curve for 

Indium at 1 g and Surrounding 
Air Pressure of 0.10 atm 
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curve. The temperature gradients are small even over the total length of the 
cylinder. Since the heater is wound around the cylinder, the temperatures are 
assumed to be symmetric about the cylindrical axia. For this reason the maxi- 
mum thermal gradient would be from the outside to the cylindrical axis, so 
that the radius is the characteristic length. The maximum thermal gradient 
along the cylindrical axis is from one end to the other end. The following 
tables gives the temperature gradients. 


Max. Thermal 
Gradient for; 

Temperature 

(°C) 

— 

Temperature 

(°C) 

AT 

(°C) 

Distance 

(cm) 

AT 

AX 

Radial 

308.353 

308.655 

0.302 

0.307 

0.98 

Longitudinal 

308.262 (In) 

308.289 (In-Pb) 

0.027 

2.35 

0.011 

Across 

In--In-Pb (ave) 
Interface 

308.506 (In) 

308.530 (In-Pb) 

0.024 

0.782 

0.031 

Across 

In--In-Pb (max) 
(next to axis) 

308.440 (In) 

308.475 (In-Pb) 

0.035 

0.782 

0.045 


1.7 CONVECTION SENSITIVITY CALCULATIONS 

This section is subdivided into three parts: (1) equations for the various 
dimensionless numbers; (2) air convection on the outside of the apparatus; and 
(3) indium and indium-lead convection inside the cylinder. 


1.7.1 Dimensionless Numbers 


The following dimensionless numbers are calculated throughout this 

section: 


Grashof number, Gr 


(Bouyancy Forces) 
(Viscous Forces) 


Prandtl number, Pr 


Momentum Transport 
Heat Transport 


16 


LOCKHEED • HUNTSVILLE RESEARCH A ENGINEERING CENTER 














LMSC-HREC TR D496846 


Rayleigh number, Ra 


(Buoyancy Forces) x (Moment Transport) 
(Viscous Forces) x (Heat Conduction) 


Nusselt number, Nu 


Convective Heat Transfer 
Conductive Heat Transfer 


Marangoni number, Ma 


Surface Tension Driving Forces 
Viscous Forces 


(Buoyancy For ces )x (Momentum Transport) 
Solutal Rayleigh - (Viscous forces) x (Mass Diffusion) 

number, Ra 


Schmidt number. Sc 


Kinematic Viscosity 
Mass Diffusivity 


Critical Rayleigh 
number, Ra^ 


Point at which convection begins to augment 
mass diffusion or thermal conduction; i.e., 
mixing begins 


The equations for each of these dimensionless numbers are: 

V 



Ra = Gr Pr 



Nu = 0.5Z5(Ra) 
do 


1/4 for convection external to a 

horizontal cylinder (Refs. 5,6 and 7) 


Ma = 


■5f 


ATH 

piTo 


Temperature gradient induced (thermal or therrnogravitational) 
Concentration gradient induced (solutal or solutogravitational) 
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where 


H 

T 

AT 

g 

P 

M 

u 

a 

a 

D 

C 

AC 

k 

(3 

h 


Ra 


P^gAC 

Di/ 


Sc 



= characteristic length* 

= temperature 

= incremental change of temperature 
= gravitational acceleration 
= density 

= dynamic viscosity 
= kinematic viscosity 
= surface tension 
= thermal diffusivity 
= mass diffusivity 
= concentration of solute 

= incremental change of solute concentration 
= thermal conductivity 
= specific heat 

= volumetric thermal expansion 
= heat transfer coefficient 


The following two-time constants are calculated for indium and indium-lead. 

* 

t , = — , momentum diffusion time (3250 sec of In Pb) 
md 1 / 

t^ , = — , thermal conduction time (25.3 sec for In R») 

td oc 

1.7.2 Air Convection 

Li order to compare the heat transfer and the convection of the air at 
various levels of gravitational acceleration, the Grashof, Rayleigh and Nuoselt 


>(< 


H = 2.54 cm. 
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numbers are calculated for various g levels. The following table (from 
Refs. 6 and 7) shows the dimensionless numbers for the air surrounding the 
cylinder under various gravity conditions. 


Acceleration 

Grashof 

Number 

Rayleigh 

Number 

Nusselt 

Number 

Ig 

4.37 X 10^ 

2.98 X 10^ 

12.26 

O.lg 

4.37 X lO'* 

2.98 X 10^ 

6.90 

lO'^g 

4.37 X 10^ 

2.98 X 10^ 

2.18 

o 

1 

4.37 X 10^ 

1 

2.98 X lO" 

1.23 


A ground test was made with 0.10 atmospheres of pressure surrounding 
the mixing cylinders and inside the container. In ordex to make a comparison 
between the actual ground experiment and the computer model, a Nusselt number 
was calculated for 0.10 atmospheres. These calculations resulted in; 


1 g and air 
at 0.10 atmo 
sphere of 
pressure 


Gr = 

4.36 X 10^ 

Nu (0. 1 atm) 

Pr = 

0.682 

Nu (1.0 atm) 

Ra = 

2.98 X 10^ 

Nu (0.01 atm) _ 

Nu = 

3.88 

Nu (1.0 atm) 


The Nusselt number was used to modify the air parameters for conductive 
and convective air under these conditions. This was done by multiplying 
the thermal conductivity by (1 + 3.88) or k' = k (1 + 3.88). 
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1.7.3 In-Pb Convection Sensitivity 

A simple, reliable estimate of whether convection exists, and its 
strength, can be obtained from dimensional analysis. This is especially true 
for thermogravitational flows which have been studied extensively (both theo- 
retically and experimentally) for the past 30 years (Ref. 8). The analysis of 
thermocapillary flows, however, is not as extensive due to the difficulty of 
suppressing the usually dominant buoyancy forces during terrestrial experi- 
ments. The following sections describe: analysis of the magnitude of thermo- 
gravitational, solutogravitational, thermocapillary and solutocapillary con- 
vection expected in rocket and ground testing; and a means for estimating 
mass transfer effects via heat and mass transfer analogies. 

• Thermogravitational Dimensional Analysis 

A schematic of the mixing experiment liquids is shown in Fig. 4. Di- 
mensional analysis of this configuration requires that values of the Rayleigh 
and Marangoni numbers be calculated for both terrestrial and rocket conditions 
and then compared to the critical Rayleigh and Marangoni numbers. The 
aspect ratio (length-to-diameter) and Prandtl number of the processes are 
also important parameters. Ignoring the Marangoni effect for the present, 
the dimensional analysis for thermal buoyancy convection is presented in this 
section. 

For heating from below or above (vertical processing) the Rayleigh num- 
ber is defined as; 
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g = 1.0 X 10"^ 

Rigid or ^ 

Free Boundaries 


Liquid 

j 

I 

1 

Liquid 

1 

In-Pb (20 wt %) 

1 

1 

In 

^C 

■ — 

I 

1 

_J 


_1 



L 


Thermal Conditions: Radial Temperature Gradient = 0.98 °C/cm 

Axial Temperature Gradient = 0.011 °c/cm 


Dimensions: L = <i,35 cm 

= 0.6l cm 


Fig, 4 - Physical Model and Coordinate System 
for InPb Experiment 


21 

LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 




) 

i 


f 


i ’■ 





whoro 




LMSC-HREC TR D496846 


('^Hot ■ '^Cold) ^ 

“■'l - wi 


P = thermal expansion coefficient 
g = gravity level 
V = kinematic viscosity 
a = thermal diffusivity 


The mean gravity expected during the rocket coast phase is approximately 
-4/2 

10 ~ cm/sec ). Thus during vertical processing, ll^L.j£a,rth ~ 

I ”*4 c 

2,7 and I Rocket “ critical Rayleigh number, depends 

on the heating direction, aspect ratio (7), and Prandtl number (Pr). For the 

mixing experiment, Pr = 8 • 10 ^ and 7 = l/D = 3.8. Thus Ra^^ = 1.6 • 10^ for 

heating from below (Ref. 9) and convection is not even likely in ground proces- 
c 

sing. Raj^ - 00 for heating from above which means that convection is theo- 
retically always absent. Actually, however, finite lateral temperature gradients 
exist (side heating or cooling); and Ra for side heating is zero (some flow 
always exists). In side heating, however, the heat and mass transport may 
not be increased until a second critical Rayleigh number, Ra^. is reached. 

For side heating (horizontal processing) the Rayleigh number is defined 

as, 

and the critical Rayleigh number for convective heat and mass transport is de- 
fined as (Ref. 10). 

Ra = 592 rir -1 = 3.8 • 10^ 
w 

2 -2 
For ground processing, Ra^ = 1.1*10 ; while for the rocket, Ra^ = i 1 • 10 
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Thus, utilizing a rocket mean gravity level of 10 g^, the Rayleigh 

numbers are: 



Ground 

Rocket 

Ra^/Rau 

2 • lO"'^ 

2 • 10'® 

Ra /Ra 
w w 

3 • 10"^ 

3 • 10'® 


Therefore, thermal buoyancy convection is unlikely during ground or micro- 
gravity processing. 

• Solutogravitational Dimensional Analysis 

For freezing from below or above (vertical processing) the solutal 
Rayleigh number is defined as: 

^ Q g AC 

= — w --- 



where 



C = concentration of solute* 


3) = molecular diffusivity of solute 


-4 

Utilizing 10 g^ for rocket processing, solutal Rayleigh numbers for vertical 
processing are: 


Earth: Ra^^ = 3.1 • 10^ 

Rocket: Ra'j^ = 3.1 • 10^ 

Q/ 

The critical solutal Rayleigh number, Ra^, depends on the freezing 
direction, aspect ratio, 7 , and Schmidt number, Sc. For the mixing experiment. 




AC = 0.2 wt % Pb, L = 1.18 cm 


23 


LOCKHEED HUNTSVILLE RESEARCH k ENGINEERING CENTER 


T 



LMSC-HREC TR IW 96846 


Si' = 170 and V= l/D = 3.8. Thus RSj^ = 5 • 10^ for freezing from above 
(In-Pb segment in upper half of cylinder) (Ref. 9) and solutal convection is 
probable in both ground and rocket processing. Ra^^ = oo for freezing from below 
(In-Pb segment in lower half of cylinder). This means that solutal convection 
is not possible in the latter configuration. . 

For side heating (horizontal processing) the folutal Rayleigh number is 
defined as; 


Ra^ = ^^g(AC) 

and the critical Rayleigh number for convective heat and mass transport is 
defined as (Ref. 10). 


Ra^ = 592 yVr^ = 3.8 ' 10^ 

For ground processing, let Ra^ = 4 * 10®; while for the rocket, let Ra^= 4 • lo'*. 
Thus, utilizing a rocket mean gravity level of 10 ^ gg, the solutal Rayleigh 
numbers are; 



Ground 

Rocket 


6.2 • lo'^ 

6.2 

Ra' /Ra 
w w 

10.5 • 10^ 

10.5 


Thus, convection is likely during freezing from above processing on the 
ground and on the rocket; while horizontal processing will also yield convective 
influences during ground and rocket processing also. 

# Thermocapillary Dimensional Analysis 

Surface tension gradients may be an important source of convection in 
melts whenever a liquid-fluid interface is present (Ref. 11). 
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As with thermogravitational convection, thermocapillary flow will be 
laminar with small surface tension driving forces, oscillatory at moderate 
driving forces and turbulent at large driving forces. The ratio of capillary 
driving force to resisting viscous force is expressed by the Marangoni 
number defined as. 




^ ATL 


for temperature gradients parallel to the liquid-fluid interface and as 


Ma 

w 


dcr 

dT 


AT w 


/>VOt 


for temperature gradients perpendicular to the liquid-fluid interface. The 
symbol cr represents the surface tension and da/dT = -0.08 dyn/cm°C for 
InPb(Ref. 12). 

A critical Marangoni number has been established for Ma (Ref. 13), but 

Q 

not for Ma, . Values of Ma. , Ma and Ma are shown below. 

L L w w 

Maj^ = 2.85 

Ma^ = 1.95, Ma^ = 70 

Even though no firm critical Ma^^ is known, values of Maj^ > 7 exhibited 
oscillatory laminar flow and Ma^^ > 140 indicated turbulent flow for molten 
silicon (Ref. 11). It is therefore evident that significant thermocapillary flow 
should not occur in the mixing experiment. 

It should be mentioned at this point, however, that little or no data exist 
on values of do/dT in the presence of dopants or of oxide films (see Appendix 
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i 


C, Ri'f. 14). Nuither the preceding Marangoni numbers nor the analysis by 
Chang and Wilcox (Ref. 11) account for these phenomena. Thick oxide films 
may completely suppress thermocapillary convection. 

• Solutocapillary Dimensional Analysis 

No explicit data exists on values of da/dC for the In-Pb system. 
Esimating da/dC from Hoar and Melford (Ref. IZ) the solutcd Marangoni 

I I 

numbers, Ma, and Ma , are as follows: 

w 

• For concentration gradients parallel to the liquid-fluid interface, 


I& 

■ Pua 


= 1.7 • 10 


• For concentration gradients perpendicular to the liquid-fluid interface. 


da 

Ma' =dC 


ACD 


w 


= 2.2 • 10 


Thus solutocapillary convection may be significant on the rocket. 

• Heat and Mass Transfer Rates 

Convection affects solidification from the melt primarily through its 
influence on (Ref. 15): (1) the amount of dopant or impurity delivered to the 
vicinity of the solid-liquid interface; i.e., the rate of mass transfer; and (2) 
the temperature fluctuations cuused near the growth interface (heat transfer). 
A convenient treasure of the convective augmentation to heat transfer is the 
Nusselt number, Nu. It represents the ratio of total heat transfer to conduc- 
tion heat transfer and is defined as 

Nu = hL/k 
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where h is an overall coefficient of heat transfer. It is related to thermal 
and concentration driving forces as follows 


C2 c 

Nu = Cj Ra + Cj Ma + C 


'^6 '^8 
- Ra° + C_ Ma® 

3 I 


where the C^'s are constants. 

The constants, C^, sometimes incorporate the effects of geometry via 
the aspect ratio, 7 , and the effects of relative diffusivity via the Pr and Sc 
numbers. The Prandtl number is the ratio of momentum diffisivity, 1 / , to 
thermal diffusivity, a . For solutal convection, the Schmidt number, Sc, rather 
than Pr is important. The Schmidt number is the ratio of momentum diffu- 
sivity to molecular diffasion {Sc =v/5)). 

Transport by convection tends to dominate as Pr or Sc becomes large. 
These two parameters are also important when estimating mass transfer rates 
from measured heat transfer rates and vice versa (Ref. 17). The degree of con- 
vective augmentation to mass transport is reflected by the Sherwood number, 

Sh, which is analogous to Nu. For small flow rates, the Sherwood, Nusselt, 
Schmidt, and Prandtl numbers are related as follows 


Sh 

Nu 



n 


where n « l/Z for laminar flow and n « l/3 for turbulent flow. The analogy 
between heat and mass transfer can be utilized to estimate mass transfer rates 
from computed heat transfer rates or vice-versa. Furthermore, the Nusselt 
and Sherwood numbers can be related as follows for vertical processing 
(Ref. 16): 
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Nu^ = 1 + .14 


Ra^ + Le 


-16/n 1 3/8 


Shj^ = 1 + (Nuj^ - 1) Le 


6/ll 


where Le = Lewis number = Sc/Pr. For freezing from above in the mixing 
experiment 


Le = 2.5 • 10^ 

= 1.04 
Sh^ = 17.6 

Thus convective mixing, generated primarily by concentration gradients, 
should lead to significant transport of solute relative to the rate of pure 
diffusion, even though convective heat transfer augmentation is minor (Nu«l). 

• Conclusions 

Sensitivity calculations indicate that thermogravitational convection will 
not be important in ground or rocket processing. Solutogravitational flow should 
be significant even for rocket processing. Thermocapillary flow should be 
negligible in both ground and rocket testing. Solutocapillary convection may also 
be significant, if the Pb-In concentration gradient were to exist at a liquid- 
vapor interface. 
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1.8 CONVECTION MODELING 

Results in this subsection are directed toward the theoretical prediction 
of the mixing experiment results. The ability to analytically model the ex- 
periment is essential for analysis of the sounding rocket flight data. The 
following is a description of the configuration analyzed and of some flight 
results that have been obtained. 

The basic portion of the experiment consists of three orthogonal cylinders 
each of which contain indium and an indium-lead alloy. The metals are melted 
and held in the molten state for about four minutes. A complete description 
of the cylinder design, heat-up etc., is given in Section 1.1, For the purpose 
of analysis the configuration in Fig. 5 is used. 

In Fig. 5: 

w = mass concentration of indium 

T - temperature 

f = jitter component of gravitational field 

g^ = earth gravity 

A = wavelength of "sawtooth" g-jitter 
The definitions of other variables and parameters in Fig. 5 are obvious. 

The gravitational field and the time variations shown are the result of 
studying the accelerometer data from the 4 October 1974 sounding rocket 
flight. The "sawtoo ..." jitter with an amplitude and frequency roughly equal 
to the average of the flight data was chosen to represent the g-field variations. 

The Lockheed Convection Analysis program was used to predict the 
mixing for the configuration in Fig. 5. The cylinder could not be analyzed for 
the situation given in Fig. 5 because the flow field is three-dimensional. Thus 
a rectangular box was analyzed to obtain some qualitative results. The phe- 
nomena assumed to be driving the mixing are diffusion and buoyancy. 
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For the gravity levels used, the analysis indicates that mixing may 
occur as far as 0.5 cm from the initial indium/indium-lead interface as a 
result of buoyancy. The mixing is expected to be approximately as shown 
in Fig. 6. 
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Section 2 

POSTFLIGHT ANALYSES 


2.1 INTRODUCTION 

This section contains results of postflight convection sensitivity- and 
modeling analyses. These studies were conducted to assist the Principal 
Investigator to interpret rocket flight results. The rocket flight results, 
both estimated g-levels and mixing, are shown in Figs. 7, 8 and 9* Results 
show that no flow was experienced in the sample mounted parallel to the 
rockets longitudinal axis. One of the samples mounted in a r>lane perpendicular 
to the rocket's longitudinal axis experiencer* a very small motion of the inter- 
face. The other sample (in the same plane) experienced flow down one side 
of the sample container. The maximum gravity levels were estimated to be 
2.10 g (Ref. 19). Thus significant mixing and complete ./ orientation of the 
lead phase occurred in Flight Sample No. 6, but no significant mixing occured 
in the other two flight samples. 

2.2 CONVECTION SENSITIVITY 

The convection sensitivity analysis of Section 1.7 has been updated to 
account for actual rocket acceleration levels and for an important, recent 
contribution (Ref. 18) to the state-of-the-art in the theory of convection onset. 
Nield (Ref. 18) has evaluated the effects of nonlinear, vertical density pro- 
files on the onset of convection (critical Rayleigh number, Ra^). In particular, 
he predicts that a step function change in density at the mid-height of the 
container will reduce Ra^ to 60% of its value for a linear density profile. 
Applying this result (in conjunction with the sensitivity criteria of Section 
1.7. 3) to the rocket results (as given by Table 3), indicates that rocket 
sample No. 6 should have indeed showed the greatest tendency toward mixing. 
Furthermore, it is important to note that the effect of step-function density 
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Sample No, 1 
’ (Dark Area is In-Pb) 






Rock et 

Longitudinal 

Axis. 

Arrow is 
I Directed 

I Toward Nose 
of Rocket . 




I 


Fig. 7 - Orientation of Sample No. 1 
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Fig. 8 - Orientations of Sanaple Nos. 2 and 6 
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c y c 

profiles on Ra is unknown. Therefore Ra /Ra »-l (which indicates 
w w w 

mixing) shown in Table 3 may not be valid. This further reinforces that 
mixing should have occurred in rocket sample No. 6 and not in the other 
rocket samples. 

2.3 CONVECTION MODELING 

The Lockheed Convection Computer program was used to model and 
analyze the rocket results. The initial and boundary conditions utilized for 
these models are identical to those of Fig. 5 except for gravity. The model 
results and g-levels used are shown in Figs. 10, 11 and 12. These figures 
show that no significant mixing occurred in any of the rocket samples. The 
melt velocities computed were less than 0.001 cm/sec for all the cases. Thus 
the extreme mixing demonstrated in Flight Sample No. 6 must have been due 
to surface tension forces or buoyancy forces larger than used in the computer 
analysis. 

The computer analysis is limited in accuracy by its inability to account 
for three-dimensional flows (two-dimensional flow in rectangles modeled). 

It is believed that gravity levels shown for Sample No. 6 in Table 3, rather 
than those shown in Fig. 12 (for Sample No. 6), are more characteristic of 
actual flight conditions and would result in model results consistent with the 
flight results shown in Fig. 9. 

Surface tension forces could also have given rise to mixing if the In/ 
In-Pb interface ever existed along a liquid/ vapor free surface. This could 
only occur if the melt pulled away from the cylinder wall at the mid-length 
of the container. This is an unlikely occurrence. 

In conclusion, it is felt that mixing in Flight Sample No. 6 was caused 
by gravity levels and directions not currently accounted for (in conjunction 
with solutally induced density gradients). Therefore, further analysis of 
flight accelerations and additional computer model simulations (using other 
gravity conditions) are recommended. 
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Single Heater Assembly (from MSFC Dwg. 95Mil309) 




NOTES: 


1. Remove all burrs and break sharp edges. 

2. Material, Ai Aly. 6061-T6-DER, WW-T-700/6. 
— 0.035 wadi thickness 
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ceramic 12 hours at 400°F. Then apply 
»xy 2651 with No. 9 catalyst and cure 
room temperature or 1 hour at 100®F. 

of Aerocoax 1 hr-040Bl6.3. Close 
ckwise on . 

51 epoxy with No. 9 catalyst; Cure for 
room temperature or 1 hour at 160 F. 

) and for fit. Installation will be 
, MSFC. 

or assembly @ will be installed by 


MSFC-SPEC-135, Class 
t, visual inspection only. 


U, sealed 
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